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Abstract A statistical study has been carried out of the relationship between 
plasma flow Doppler velocities and magnetic field parameters during the emer- 
gence of active regions at the solar photospheric level with data acquired by the 
Michelson Doppler Imager (MDI) onboard the Solar and Heliospheric Observa- 
tory (SOHO). We have investigated 224 emerging active regions with different 
spatial scales and positions on the solar disc. The following relationships for the 
first hours of the emergence of active regions have been analysed: i) of peak 
negative Doppler velocities with the position of the emerging active regions on 
the solar disc; ii) of peak plasma upflow and downflow Doppler velocities with 
the magnetic flux growth rate and magnetic field strength for the active regions 
emerging near the solar disc centre (the vertical component of plasma flows) ; iii) 
of peak positive and negative Doppler velocities with the magnetic flux growth 
rate and magnetic field strength for the active regions emerging near the limb 
(the horizontal component of plasma flows) ; iv) of the magnetic flux growth rate 
with the density of emerging magnetic flux; v) of the Doppler velocities and 
magnetic field parameters for the first hours of the appearance of active regions 
with the total unsigned magnetic fiux at the maximum of their development. 
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1. Introduction 



Solar magnetic flux emerges at various spatial scales (e.g. Parnell et al. 2009). 



Golovko ( 1998 ) revealed a power-law relation between the maximum magnetic 
fliDc and the lifetime of active regions. Otsuji et al. ( |2011 ) found a relation 
between the flux growth rate and the maximum magnetic flux in the form of a 
power-law function. 
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Measurements of Doppler velocities indicated that there is a vertical plasma 
upflow up to 1 km at the polarity inversion line when active regions emerge in 



the solar photosphere ( 


Brants, 1985a', 'Brants, 1985b 


Tarbell et al., 1989 


Lites, 


Skumanich, and Martinez Pillet, 1998;,Strous and Zwaan, 1999| Kubo, Sh 


imizu. 


and Lites, 2003 Guglielmino et al, 2006 Grigor'ev, Ermakova, and Khlystova, 


2007 Grigor'ev, Ermakova, and Khlystova, 2009 


I. 


Grigor'ev, Ermakova, and 


Khlystova (2007) revealed high Doppler velocities of about 1.7 kms ^ at the 



beginning of the emergence of the powerful active region NOAA 10488 at heli- 
ographic coordinates NOB E31 (Bo+4.9°). The observed Doppler velocities have 
not been compared with the parameters of the magnetic fields. 

The cold and dense plasma contained inside magnetic flux emerging into the 



2 kms ^ at the photospheric level (Gopasyuk, 1967 Gopasyuk, 1969 Kawaguchi 


and Kitai, 1976; Bachmann, 1978 Zwaan, Brants, and Cram, 1985 Brants, 


1985a 


|Brants, 1985b, ,Brants and Steenbeek, 1985| Lites, Skumanich, and Mar- 


tinez 


Pillet, 19981 Solanki et al, 2003 |Lagg et al, 2007 |Xu, Lagg, and Solanki, 


20101 


. Observations of pores revealed a statistical relation between the Doppler 



velocity of plasma downflow and the magnetic field strength. |Bonaccini et al. 
(1991), studying a large pore, found that the strong magnetic fields correlate 
linearly with downflow D oppler veloc ities in the form B ~ 500 V, where B is in 
Gauss and V is in kms~^. Keil et al. ( jl999, ) found that downflow Doppler veloc- 
i ties in side pores correlate positively with the magnetic field strength. |Cho et aJ] 
( |2010 ) analysed small pores less than 2", which were not related to magnetic 
flux emergence, by using high spatial resolution data from Hinode. The authors 
revealed a negative correlation between positive Doppler velocities and the mag- 
netic field strength, i. e. small downflow Doppler velocities corresponded to strong 
magnetic fields. Grigor'ev, Ermakova, and Khlystova (2011) have shown that in 



the emerging active region during the formation of pores the plasma downflow 
Doppler velocity increases linearly with the magnetic field strength. 

Estimations of horizontal velocities in the emerging active regions were per- 
formed indirectly. A wide range of velocity values from 0.1 to 5.5 kms~^ were 



obtained by analysing displacements of magnetic field structures (Frazier, 1972 



Schoolman, 1973| [Harvey and Martin, 197"3l |Chou and Wang, 1987| |Barth anc 



Livi, 1990 Strous and Zwaan, 1999 Hagenaar, 2001 Grigor'ev, Ermakova, and 



Khlystova, 2009 Otsuji et al., 2011). Considering the separation of opposite 



polarity poles in 24 bipolar pairs Chou and Wang (19871 found no relation 
of horizontal velocities neither with the average magnetic field strength nor 
with the total unsigned magnetic flux. Otsuji et al. (2011) studied 101 emerging 



flux regions with various spatial scales. They showed a power-law relation with 
negative index between the horizontal velocities of the separation of opposite 
polarity poles and the maximum value of magnetic flux. 



Khlystova (2011) presented a statistical study the photospheric Doppler ve- 
locities in 83 active regions with a total unsigned magnetic flux above 10^^ Mx. A 
centre to limb dependence of negative Doppler velocities was established, which 
indicates that the horizontal velocities of plasma outflows are higher than the 
vertical upflow velocities during the flrst hours of the emergence of active regions. 

In this article we study the statistical relations between Doppler velocities 
and magnetic field parameters during the emergence of active regions. 
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2. Data Analysis 



We used full-disc SOHO/MDI data (Scherrer et al, 1995). The temporal reso- 



lution of the photospheric magnetograms and Dopplergrams is 1 minute, while 
that of the continuum is 96 minutes. The spatial resolution of the data is 4". 



Magnetograms with a 1.8.2 calibration level were used (Ulrich et al., 2009). 

We have cropped a region of emerging magnetic flux from a time sequence 
of data taking into account its displacement caused by solar rotation. This 
displacement was calculated by the differential rotation law for photospheric 



magnetic fields (Snodgrass, 1983) and the cross-correlation analysis of two mag- 
netograms adjacent in time. A precise spatial superposition of the data was 
achieved by cropping fragments with identical coordinates from simultaneously 
acquired magnetograms and Dopplergrams. For correct processing of the data, 
we chose the cropped region in a way that it excluded the area outside the 
limb. For the Dopplergrams the contribution of solar differential rotation and 
other factors distorting the Doppler velocity signal were removed by using the 



technique described in Grigor'ev, Ermakova, and Khlystova ( 2007 ) 



The sequence of the fragments obtained has been used to calculate the tempo- 
ral variations of the parameters under study. The calculation area was limited to 
the region of the emerging magnetic flux. The boundary of emerging flux region 
was visually inspected. 

The total unsigned magnetic flux of active regions was calculated inside iso- 
lines ±60 G taking into account the projection effect and assuming that the 
magnetic field vector is perpendicular to the solar surface: 

'^=sof:^, (1) 

^-^ COSdi 
i—1 

where $ is the total unsigned magnetic flux in Mx, Sq is the area of the solar 
surface of pixel in the centre of solar disc in crri^, N is the number of pixels 
where \Bi\ > 60 G, Bi is the line of sight magnetic field strength of ith pixel 
in G, 6i is the heliocentric angle of ith pixel. The maximum total unsigned 
magnetic flux ^max was determined by the maximum inflection point in the 
increase of the magnetic flux curve or by the last value of the datasets (for active 
regions passing beyond the west limb or with insufficiently downloaded data). 
The signal background that existed before the emergence of active regions was 
subtracted from the maximum total unsigned magnetic flux. Magnetic saturation 
in measurements of the magnetic field strength in SOHO/MDI was controlled 



(Liu, Norton, and Scherrer, 2007); it was attained only in three of the active 
regions considered. 

The total unsigned magnetic flux growth rate d^/dt and absolute value of the 
maximum magnetic field strength H^ax were calculated in the first 12 hours of 
the emergence of the magnetic fiux for large and small active regions or from 
beginning to maximum of the total unsigned magnetic fiux for ephemeral active 
regions. 

Photospheric plasma flows are characterised by the peak (absolute maximum) 
values of negative and positive Doppler velocities Vmax- and V„iax+ which were 
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determined in the first 12 hours of the magnetic flux emergence for large and 
small active regions or from beginning to maximum of the total unsigned mag- 
netic flux for ephemeral active regions. The magnetic flux emergence begins with 
the appearance of the loop apex where the magnetic field is horizontal; therefore, 
the time interval we analysed started 30 minutes before the appearance of the 
line of sight magnetic fields. Continuum images were used to control the possible 
formation of sunspots with Evershed flows in the first hours of the emergence of 
active regions to eliminate their contribution. The maximum Doppler velocity of 
these flows in SOHO/MDI data with low spatial resolution can reach 2 kms^^ 
(Bai, Scherrer, and Bogart, 1998). 

The position of the active region on the solar disc is expressed by the helio- 
centric angle 9: 

6* = arcsin(r/i?), (2) 

where r is the distance from the disc centre to the place of the active region 
emergence and R is the solar radius. is also approximately the angle between 
the normal to the surface and the line of sight component of the emerging 
magnetic flux. 

We performed a regression analysis of the data. There is at least one inflection 
point in the considered dependencies. Therefore, linear or quadratic polynomials 
were used as approximation functions. The choice of the polynomial degree 
was based on the minimum values of the standard deviations. The results are 
given in the form of polynomial equations which, according to the F-statistic, 
are significant. The confidence intervals for the means were calculated with a 
confidence probability of 99%. 



3. The Investigated Active Regions 

We have studied 224 active regions which emerged on the visible side of the solar 
disc from 1999 to 2008 (see Figure [I]). The selected active regions have different 
spatial scales and emerge at different distances from the solar disc centre. They 
are isolated from high concentrations of existing magnetic fields (the presence 
of single poles with total unsigned magnetic fiux not higher than 0.5 x 10^^ Mx 
in the region of the direct emergence of the active region for the first hours 
was allowed). Sufficiently complete data series with a temporal resolution of 
1 minute for the first hours of magnetic fiux emergence are available for the 
selected objects. 

The statistical analysis includes the following. 

1. According to the maximum value of the total unsigned magnetic flux (in the 
two polarities): 

• 68 ephemeral active regions (2.6 x 10^^ Mx < ^rnax < 4.07 x 10^° Mx); 

• 130 smaU active regions (4.07 x 10^° Mx < <^„iax < 8 x lO^^ Mx); 

• 26 large active regions {^max > 8 x 10^^ Mx). 
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Figure 1. Maximum total unsigned magnetic flux of the (a) ephemeral, (b) small and (c) 
largo active regions versus the heliocentric angle that corresponds to the beginning of the 
emergence of magnetic fluxes. The total magnetic flux maximum was determined by the 
maximum inflection point in the increase of the magnetic flux curve (marked by points) or by 
the last measurement (marked by triangles). Active regions with magnetic saturation are also 
marked by triangles. In active regions marked by triangles, the maximum of the total unsigned 
magnetic flux may be higher. 



The limits of the total unsigned magnetic flux for ephemeral active regions 
were taken from Hagenaar (20011. The limit between the large and small 



active regions is considered to be the maximum magnetic flux in one polarity; 
it is taken equal to 5 x 10^^ Mx (Garcia de La Rosa, 1984 Zwaan, 1987). In 



the present investigation the limit between the large and small active regions 
was chosen at the total unsigned magnetic flux in two polarities, equal to 
8 X 10^^ Mx, since all active regions with the total unsigned magnetic flux 
above this level contain sunspots and their spatial size exceeds the size of a 
supergranule. 

2. According to the distances from the solar disc centre: 

• 72 active regions near the solar disc centre {0 < 25°); 

• 98 active regions at medium distance from disc centre (25° < 9 < 50°); 

• 54 active regions near the limb {9 > 50°). 

The maximum Doppler velocities of the convection flow in the quiet Sun 
were calculated at different distances from the disc centre using SOHO/MDI 
data (Figure [2]). Dopplergrams with 1 minute temporal resolution for the period 
18-23 April 2001 were used. The peak values of negative and positive Doppler 
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Figure 2. Peak (a) negative and (b) positive Doppler velocities of the convective flows in the 
quiet Sun versus the heliocentric angle for SOHO/MDI data. 



velocities were measured within a region of size 40"x40" whose displacement 
was tracked based on the differential rotation law of Doppler structures in the 
photosphere ( Snodgrass and Ulrich, 1990 ) . The selected region was going through 
the WOO - W75 range of longitudes along the equator during the period under 
study. The contribution of solar differential rotation and other factors distorting 
the Doppler velocity signal were removed by using the technique described by 



Grigor'ev, Ermakova, and Khlystova (2007). 

In Figure [2| one observes that Doppler velocities of convection flows increases 
with the heliocentric angle. For the points with 9 > 60° there is a tendency for 
Doppler velocities to decrease. This may be due to the Doppler velocities being 
measured in higher atmospheric layers near the limb, due to the increase in the 
optical thickness, where the atmosphere is convectively stable and convection 
flows are damped. 



4. Results 

4.1. The Centre to Limb Dependence of Negative Doppler Velocities 

The peak values of negative Doppler velocities (plasma motion toward the ob- 
server) in the first hours of magnetic flux emergence of ephemeral, small and 
large active regions show an increase with the heliocentric angle (Figure |3]). The 
regression analysis of the data yielded Equations ([3| for ephemeral, Q for small, 
and ^ for large active regions: 

Kiaa;- = -211.32 - 24.7761 + 0.166*2, (3) 
Vmax- = -134.93 - 33.050 + 0.266'^ (4) 

Vmax- = -96.35 - 34.7561 + 0.226*2, (5) 

where Vmax- is the peak negative Doppler velocity in ms^^ and 9 is the helio- 
centric angle in degrees that corresponds to the position of the region at Vmax- ■ 
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Figure 3. Peak negative Doppler velocities (plasma motion towards the observer) in the 
first hours of the emergence of (a) ephemeral, (b) small and (c) large active regions versus 
the heliocentric angle that corresponds to the position of the region at this time. The solid 
line corresponds to Equations ijsjl in plot (a), ^ in plot (b) and (jsj in plot (c); dotted lines 
correspond to 99% confidence intervals for the means. 



Correlation ratios show a strong relation between the considered parameters: 
—0.85 for ephemeral, —0.79 for small and —0.79 for large active regions. The 
dependencies obtained indicate that during the emergence of active regions the 
horizontal velocities of plasma outflows {9 > 50°) exceed the vertical upflow ve- 
locities {9 < 25° ) . Generally, the Doppler velocities accompanying the emergence 
of active regions are higher than those of convective flows, however, sometimes 
they are comparable. For small and large active regions the horizontal Doppler 
velocity component substantially exceeds horizontal Doppler velocities of the 
convection flows in the quiet Sun, which are less than 1200 ms~^ (Figure [sjb, c, 
and Figure |2]a the points with 9 > 50°). 

4.2. Vertical Doppler Velocities 

For the active regions emerging near the solar disc centre {9 < 25°), peak 
negative Doppler velocities (plasma upflow) have a high dispersion and cor- 
relate neither with the magnetic flux growth rate (the correlation ratio is —0.30) 
nor with the maximum magnetic field strength (the correlation ratio is —0.22) 
(Figure |4]) . The dependencies obtained imply that the upflow Doppler velocities 
do not characterise the power of the emerging magnetic flux. 
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Figure 4. Peak negative Doppler velocities (plasma motion toward the observer) versus (a) 
the magnetic flux growth rate and (b) the maximum magnetic field strength in the first hours 
of the emergence of active regions near the solar disc centre (9 < 25°). Black symbols are for 
ephemeral active regions, red symbols for small active regions and blue symbols for large active 



regions. 
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Figure 5. Peak positive Doppler velocities (plasma motion away from the observer) versus 
(a) the magnetic flux growth rate and (b) the maximum magnetic field strength in the first 
hours of the emergence of active regions near the solar disc centre {9 < 25°). The solid line 
corresponds to Equations (jojl in plot (a) and (iTjl in plot (b); dotted lines correspond to 99% 
confidence intervals for the means. The symbolcolours are the same as in Figure [4] 



For the active regions emerging in the central part of the solar disc {0 < 25°), 
the peak positive Doppler velocities (plasma downflow) are related quadratically 
to the magnetic flux growth rate (the correlation ratio is 0.75) and linearly with 
the maximum magnetic field strength (the correlation ratio is 0.75) (Figure [5]). 
The linear dependence shown on Figure [5lb c onfirms the results obtained pre- 
viously by Bonaccini et al. (19911, Keil et al] ( |1999 ), Grigor'ev, Ermakova, and 
Khlystova (2011). The Doppler velocity values substantially exceed convection 
Doppler velocities in the quiet Sun (Figure [5] and Figure [2]b the points with 
9 < 25°). 

Regression analysis of the data yielded the following equations: 



V„,ax+ = 455.82 + 6.98 x IQ-'^'d^/dt - 1.75 x 10-^»(d$/dt) 



(6) 



217 M + 0.43 X H„ 



(7) 
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Figure 6. The peak positive and negative Doppler velocities versus (a) the magnetic fiux 
growth rate and (b) maximum magnetic field strength in the first hours of the emergence of 
active regions near the li mb {Q > 50°) . The solid line corresponds to Equations (jsl and (j9l in 
plot (a) and Equations | |10| and in plot (b), dotted lines correspond to 99% confidence 
intervals for the means. I'he symbol colours are the same as in Figure [4] 



where Vmax+ is the peak positive Doppler velocity in ms^^, d^/dt is the 
magnetic flux growth rate in Mxh^^ and Hmax is the maximum magnetic field 
strength in G. 

The draining of the cold plasma, being carried out by the emerging magnetic 
flux from the convective zone into the solar atmosphere, is considered to be the 
main reason for the plasma observed downflows. 

4.3. Horizontal Dopper Velocity 



Khlystova (2012) has carried out a detailed analysis of photospheric flows in 
four active regions emerging near the solar limb. It has been found that ex- 
tended regions of Doppler velocities with different signs are formed in the first 
hours of the magnetic flux emergence in the horizontal velocity field. The ob- 
served Doppler velocities substantially exceed the separation velocities of the 
photospheric magnetic fiux outer boundaries and most likely are caused by the 
significant component of the velocities of plasma downflow being carried out into 
the solar atmosphere by emerging magnetic flux. 

One observes in Figure [3] that the horizontal Doppler velocity component 
{9 > 50°) in the region of emerging magnetic fluxes has a wide range of values. 
For the active regions emerging near the limb {9 > 50°), the peak positive 
and negative Doppler velocities of photospheric plasma flows show a quadratic 
relation with the magnetic flux growth rate (correlation ratios 0.78 and —0.73) 
and a linear relation with the maximum magnetic fleld strength (correlation 
ratios 0.58 and —0.56) (Figure [6|. Strong magnetic flelds are oriented mainly 
perpendicular to the surface after the emergence. Therefore, the magnetic field 
strength in active regions near the limb is 2-3 times smaller than that one in 
active regions near the solar disc centre, because of the projection effects of the 
magnetic field vector to the line of sight (Figure |6]b and Figure [5]b). 

Regression analysis of the data yielded the following equations: 

Vmax+ ^ 782.27 -I- 6.68 x lO^^^d^/dt - 9.65 x IQ-^^ {d<i> / dtf , (8) 
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Figure 7. Magnetic flux growth rate versus tiie magnetic field strength in the first hours of 
the emergenc e of active regions near the solar disc centre {6 < 25°). The solid line corresponds 
to Equation | |12| l; dotted lines correspond to 99% confidence intervals for the means. The 
symbol colours are the same as in Figure |4] 



V„,ax^ = -890.85 - 5.95 x lO'^^d^/dt + 9.94 x 10"^^(d$/rfi) 



(9) 



= 630.53 + 0.85 x H„ 



(10) 



F„a.- = -770.35 - 0.71 x H^ax, 



(11) 



where Vmax+ V^ax- ai'Sj respectively, the peak values of positive and neg- 
ative Doppler velocities in ms"^, d^/dt is the magnetic flux growth rate in 
Mxh^^ and Hmax is the maximum magnetic field strength in G. 

One sees in Figure [5]a and Figure [6]a that the horizontal Doppler velocities 
exceed the vertical downflow Doppler velocities by a factor of 1.5-2 in the 
active regions with the same magnetic flux growth rate. This supports theoretical 
models in which the velocities of draining plasma have a significant horizontal 



component at the beginning of the magnetic flux emergence (see Shibata et at. 
(|1990D, lArchontis et al.\ (|2004|), iToriumi and Yokoyamal (|2010|), iToriumi et al. 



(2011) and others). 



4.4. The Relation between the Magnetic Flux Growth Rate and the Maximum 
Strength of Magnetic Fields 

A quadratic relation with relatively low dispersion and a high correlation ratio 
of 0.88 between the magnetic flux growth rate and the density of the magnetic 
fields for the active regions emerging near the solar disc centre {9 < 25°) has 
been found (Figure [7]). The equation derived by the regression analysis has the 
following form: 



= 3.52 X 10^-* - 7.27 x 10^^ H,nax + 7.49 x 10"H, 



(12) 



where d$ /dt is the magnetic flux growth rate in the first hours of the emergence 
in Mxh^^ and Hmax is the maximum magnetic held strength in G. 
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Figure 8. The relation between (a) peak negative Doppler velocities, (b) peak positive 
Dopplcr velocities and (c) the magnetic flux growth rate in the first hours of the emergence of 
active regions with total unsigned magnetic flux at the maximum of their development ^max ■ 
The symbol colours are the same as in Figure |4] The ^max was determined by the maximum 
inflection point in the increase of the flux curve(marked by points) or by the last measurement 
(marked by triangles). In active regions marked by triangles, the ^max may be higher. 



The dependence is theoretically expectable. The magnetic flux emerge due 
to the action of the magnetic buoyant force which is proportional to the square 
of the magnetic field density (Parker 19551. The higher the magnetic field 



strength, the higher the magnetic pressure and the lower the gas pressure inside 
the magnetic structure. The low gas pressure leads to a large magnetic buoyancy 
force and, thus, to an increase in the magnetic flux growth rate. 



4.5. The Relation of Parameters between the Beginning and Maximum 
Development of Active Regions 



Active regions emerging near the solar disc centre {9 < 25°) are considered. Peak 
negative Doppler velocities in the first hours of the magnetic flux emergence 
do not show any relation with the total unsigned magnetic flux at maximum 
development of the active regions (Figure [sja) . Comparison of the peak positive 
Doppler velocities and the magnetic flux growth rate with the maximum value 
of total unsigned magnetic flux reveals a power-law relation with high dispersion 
(Figure [sjb, c). A power-law relation, but with less dispersion, between the flux 
growth rate and the the maximum value of magnetic flux has been found earlier 



by Otsuji et al. (2011 ). In the present investigation the weak dependence between 



SOLA: khlystova.tex; 18 December 2012; 3:31; p. 11 



A. Klilvstova 



the parameters under consideration can be explained by the fact that active 
region magnetic flux emerges in separate fragments on time scales from scvciral 
hours to 5-7 days. Therefore, the parameters of magnetic fields and Doppler 
velocities in the first hours of emergence do not characterise an active region as 
a whole. 

5. Conclusions 

The study of plasma flow Doppler velocities and magnetic field parameters 
during the emergence of active regions at the photospheric level has revealed 
the following statistical relationships: 

1. Peak values of negative Doppler velocities accompanying the magnetic flux 

emergence in the first hours show an increase with the heliocentric angle 9. 
The observed relation shows that horizontal velocities of plasma outflows are 
higher than the vertical upflow velocities. 

2. The upflow plasma Doppler velocities have no relation with parameters of 
magnetic fields and, thus, do not characterise the power of the emerging 
magnetic flux. 

3. The downflow plasma Doppler velocities have high values, which substantially 
exceed Doppler velocities of convection flows in the quiet Sun. They are 
related quadratically with the magnetic flux growth rate and linearly with 
the maximum magnetic field strength. The main reason for the observed 
downfiows is considered to be a draining of the cold plasma, being carried 
out into the solar atmosphere by the emerging magnetic flux. 

4. The horizontal plasma outflow Doppler velocities arc connected quadratically 
with the magnetic flux growth rate and linearly with the maximum magnetic 
field strength. Possibly the horizontal flows are connected with draining of 
plasma which has a signiflcant horizontal component at the beginning of the 
magnetic flux emergence. 

5. There is a quadratic dependence with relatively low dispersion and a high 
correlation ratio of 0.88 between the magnetic flux growth rate and the maxi- 
mum strength of magnetic flelds. The observed relation is consistent with the 
magnetic buoyancy law. 

6. There is a weak power-law dependence of peak positive Doppler velocities and 
the magnetic flux growth rate in the flrst hours of the emergence of active 
regions with the maximum value of the total unsigned magnetic flux. 
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